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Abstract 

This is a technical document that provides supporting information and details of the 
publicly available code used for the preparation of the analysis for preprint "A search 
for heavy Kaluza-Klein electroweak gauge bosons at the LHC^ (submitted to JHEP). 

The Moses C++ framework is a project written for probing and developing new 
models for High Energy Physics processes which allows complete events to be sim- 
ulated by interface with the standard simulation program PythiaS Q. This paper 
demonstrates the usage of MoSES in a study of the nature of Kaluza-Klein (KK) 
excitations in a specific model where the SU{2) x U{1) gauge fields can exist in a 
single Extra Dimension (ED) compactified on a S^/Z2 orbifold, while the matter 
fermions and SU (3) gauge fields are localized in the 3d-brane. Using this framework, 
the events have been fully simulated at hadron level including initial and final state 
radiation. The study of particle decays was used to develop a method to distinguish 
between this Kaluza-Klein model and processes with similar final states. As a con- 
sequence, the possibility of observing and identifying a signal of the first excited KK 
state of the 'y/Z^ bosons in the LHC is also discussed. 
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1 About the project 

The Moses framework was developed during an MCnet 4-months project at University 
College London. Subsequently further development took place in Tel Aviv university. 
This paper discusses the first proof-of-concept version. The code has been developed in 
C++ using standard Gnu development tools running on Scientific Linux. It requires the 
packages ROOT LHAPDF [|], HepPDT 0] and PythiaS to be accessible. This paper 
contains two parts: (a) the development of the validation processes, and (b) the specific 
implementation of the KK process and the corresponding analysis of the events at the 
generator level. 
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2 Introduction 



2.1 Heavy Kaluza-Klein gauge bosons search at the LHC 

The LHC is expected to be ready for colliding beams at the end of 2009. It is designed to 
collide proton beams at 14 TeV, the highest CM energy ever reached in a laboratory. It 
will greatly enlarge the kinematic region for the search for physics phenomena Beyond the 
Standard Model (BSM). Several BSM theories predict the existence of other dimensions 
in addition to the usual three spatial and one time dimension. These models allow various 
particles to propagate into the extra-dimensional bulk. The TeV~^ ED model considered 
here allows the KK modes of SU{2) x U{1) gauge fields to propagate into the extra- 
dimensional bulk while restricting all the matter fermions and the SU{3) gauge fields to 
be localized in the usual 3d brane One objective of the LHC program in the context 

of these models is to search for a signal of the first excited KK mode of the SU{2) x U{1) 
gauge fields, denoted by 7* and Z*. 

A popular model that does not involve extra dimensions is also considered for compari- 
son. An extra heavy boson arising from the breaking of the Eq group is assumed. The 
signal of this extra heavy boson, denoted by Z', can demonstrate similar characteristics to 
the KK signal for which a technique to distinguish between the two models is required. 

Another popular possibility that is not considered in this paper is the RS model 
that predicts the existence of Kaluza-Klein spin-2 gravitons, denoted by G. In some models, 
the mass of first KK graviton excitation can be the same as the mass of the spin-1 '-/* /Z* 
or Z', and therefore, there is a need to identify the G signature and especially its spin. 

In this paper, the signature of the first two cases is studied where the produced bosonic 
candidates from either the Standard Model (SM) '~f/Z^ bosons, the heavy KK 'j* /Z* bosons, 
or the extra bosons Z', decay all into charged leptons. 

This KK model is particularly interesting because of the strong destructive interference 
that manifests itself at much lower invariant masses comparing to the resonance itself. 
In this model this will always occur for masses around half of the resonance mass. For 
example, with a resonance at 4 TeV, this will occur around 2 TeV. Therefore, at the LHC, 
the suppression of the cross section can be observed much earlier than the resonance itself. 
This will not happen for the various Z' possible signals that go along with the SM Z-line 
shape up to masses near the Z' resonance |TI[ . 



Powerful methods to quantify the sensitivity for BSM physics at the ~1 TeV scale are 
presented. For an observed resonance above a small SM background around ~4 TeV, a 
measurement of the charged lepton kinematic distributions can provide the discrimination 
between the two spin-1 BSM models presented in this paper. The sensitivity for this 
discrimination depends on the LHC luminosity and the collisions energy, and it should be 
possible already with an integrated luminosity of £=100 fb^^ by applying the Kolmogorov 
test on the invariant mass and the angular distributions^. This measurement can also 
quantify some interesting physical aspects of the new models and support the discrimination 
between them to some extent. The main observable that is extracted, in this context, is the 
forward backward asymmetry^ of the angular distributions. Although it is not shown here, 
this measurement can enable the classification of the resonance as spin-1 or, for instance, 
spin-2 to some extent^. Further consistent experimental studies on Z\ 'y* /Z* and the G 
possible signals in the LHC can be found elsewhere |T2|-pr5[]. 

2.2 Project overview 

Many processes are implemented in standard event generators, however, there are numerous 
specific BSM processes which are not modeled. This framework can be used to integrate 
various new BSM processes with the standard event generator PythiaS. New processes 
can of course be directly implemented in PythiaS where there is a special interface class 
to that job. In some cases it is natural and convenient to integrate the new model. In 
other cases, it may turn out to be more complex. In that light, the MoSES framework 
has some advantages. It allows, for instance, to utilize the interfaces with the standard 
HepPDT, LHAPDF and ROOT tools. In addition, the MoSES structure enables to form 
new independent modules that assemble the new model, as it was originally built to allow 
programming of complicated models with several non-standard aspects like the KK model. 
In the new BSM models, we include new models where two initial state particles interact 
and create a maximum of three final state particles, ze, 2 — 1, 2 — 2 and 2 — )■ 3 hard- 
processes described by a differential cross-section function, as these are the constraints 
introduced by PythiaS. Within MoSES, these processes (cross-section functions) can be 

•^The Kolmogorov test is applied between the pseudo-data samples and the MC reference samples. 
^There is a big difference between the KK and Z' asymmetry where the KK behave similar to the SM. 
^For that purpose, the full angular distribution, of both decay angles, can be used fl^ . 
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both analyzed independently and interfaced with PythiaS to generate the corresponding 
fully simulated physics events. 

PythiaS, the new generation of the commonly used event generators, Pythia, is a pow- 
erful and convenient tool. Of the most important feature in this context, is the care that 
has been taken to allow user supplied hard subprocesses to be quickly and fully integrated 
into the PythiaS framework. This was the reason the MoSES framework was developed on 
top of PythiaS, and will be enhanced with more interfaces and examples in the future. 

It is suggested that introducing new processes will be forehanded by a validation pro- 
cedure using analogous processes. Since the presented KK implementation has many com- 
mon characteristics with the corresponding SM implementation, a validation procedure 
was adopted. The internal PythiaS SM process was reproduced as a user process but while 
adjusting the helicity-amplitude formalism intended for the KK scenario. The helicity am- 
plitude was sent to PythiaS using its standard SigmaProcess interface class. The 2 — 2 
interface was used although it is also possible to utilize the 2 — )• 1 production and sub- 
sequent 1 — )■ 2 decay. The externally generated events were then compared with those 
from the self internal implementation subprocesses of PythiaS. Throughout the following 
sections, the "external / internal processes" nomenclature will be used in this context. 



Throughout the presented work the MRST20011o parton distribution set was used. 



3 Validation processes 

In the following section two methodical case studies are presented to validate the im- 
plementation of the external processes against PythiaS's internal processes; (a) the SM 
e^e~ —7- 7 —7- fi'^fi^ at low energies where the contribution of the boson is negligible 
and (b) the SM qq ^ /Z^ ^ I'^l" at higher energies where the contribution of the Z^ 
is dominant. Apart from the validation objective of these two examples, it is worthwhile 
going through these cases in some detail since it will also serve to rigorously explain the 
formalism used in the remainder of this paper. 

3.1 The SM s-channel photon exchange e^e~ — )- 7 — )- fjL'^fjL~ 

In this case the PythiaS output of the estimated total cross section values using the external 
implementation and using PythiaS's internal scheme are compared. The integrated total 
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cross section <J{S), 

.W = ^^^. (1) 

3 s 

where h = c = 1, can be realized by averaging the differential cross section over all the 
incoming helicity states, (corresponding to an unpolarized beam), by summing over all the 
outgoing helicity states, (corresponding to an unmeasured final polarization state) and by 
integrating over the solid angle dQ. The differential cross section itself in terms of the 
helicity of the incoming electron Ag- and the helicity of the outgoing muon A^- is 

da {s, cos 9) \p..p.. 2 



dn As (25e- + 1) (25e. . _ 1 _ 1 



7^ E E (l+4Ae-A,-cos^)^ (2) 



where, y/s is the collision's CM energy, the quantities Cg are the charges (in units of the 
proton charge) of the leptons, see the corresponding tree-level diagram in Fig 0. 




Figure 1: The s-channel photon exchange Feynman diagram 



The number of possible incoming helicity states {23^+ + 1) (2S'e- + 1) is expressed here 
in terms of the spins of the colliding particles. From helicity conservation, it is sufficient to 
sum only over the helicity states of the incoming electron and the outgoing muon. Finally, 
the polar angle 6 is the angle of the outgoing /i~ relative to the incoming e~ (with the 
azimuthal angle 0, distributed uniformly). If one is to use the differential cross section 
from Eq ^ in the SigmaProcess 2 — )■ 2 interface class of PythiaS, it is expected to be 
given in the Mandelstam variable t related to cos 6 as 




The transformation introduces an extra | factor in the cross section. Note that by doing 
so, the differential cross section is now dimensionally different than the previous by ^ and 
therefore has the units of GeV~^, as required by PythiaS, 

da 2 da , , , 

1^ = -71 a- 

dt s d cos & 
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After integrating over the azimuthal angle 0, the correct function given to PythiaS is 

Ag dig- Ay^ 

One can use Eq|^ to generate events with low CM (CM) energy - at 20 GeV, far enough 
below the pole. In this case it acceptable to ignore at the first approximation the 
contribution. All parton-level switches are turned off so the events are generated up to the 
level of the hard-process {ie, no parton showering, hadronization, kinematic cuts, etc.) 

For comparison, the generation of the hard-subprocess is performed using the above but 
also repeated with the same run conditions and with the same sample size of IM events 
but simply calling the internal process. In fact, in PythiaS this is usually performed in 
two steps; the 2 — ?■ 1 production and then the 1 — 2 decay. However, the statistical error 
involved in the 2—7-1 production is very small since the angle is already integrated out 
and therefore, all phase space points will be the same and evaluate to the same value. 
Therefore, the error which is subjected to roundoff errors, is unrealistic small and can be 
ignored. For this reason, there's also a less familiar 2-^2 implementation (in PythiaS) 
which was chosen for comparison since in this way one can get different cross sections 
depending on the angle selected on event by event basis and that leads to different event 
weights, which implies a realistic statistical error but also a lower efficiency. The 2 — )■ 2 
PythiaS internal implementation had to be slightly modified to fix the final state at fJ^^fi'- 

The external function (see Eq ^) describes the same process and was implemented as a 
2—7-2 with the intermediate photon included explicitly in it. The results are summarized 
in Table |l| where the agreement between these two values is to within ~l-sigma. 

Table 1: Cross section statistics as given by PythiaS for the two IM samples (internal & 
external) of e+e" — )■ 7 — > /i events at 20 GeV CM energy. 



Source code 


Estimated o"totai[nb] 


Internal 


0.2391 ±0.0001 


External 


0.23S9± 0.0001 
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3.2 The SM s-channel 7/Z° boson exchange qq j/Z^ — )- /+/ 

As with the previous example, this process can also provide useful information when vali- 
dating our mechanism. As before, all parton-level switches are turned off. The first stage 
is to compare the estimated total cross section values of the external implementation with 
the internal implementation from PythiaS for the hard process products only. The second 
stage is the detailed comparison of the differential cross sections. 

Since initially the focus was only on low masses, only the 7 contribution was taken into 
account. At higher energies we must fully take into consideration the contributions from 
the Z° boson and 7 — Z° interference terms. By looking only at the leptonic final states 
which are produced by the s-channel processes we can ignore the more general qq — )■ // 
processes produced by both the s- and the t-channel exchange where / can be any fermion. 
This is experimentally useful, since the background from the SM QCD interactions were 
quark and gluon final states to be considered would be very high. 

3.2.1 Hadronic level: ^ 7 — )■ l+l~X 

The Drell-Yan cross section for qq annihilation to a charged lepton pair via an intermediate 
massive photon can be easily obtained from the fundamental e~^e~ — 7 — fi^fi" cross 
section by the introduction of the appropriate color factors and by replacing the electron 
charge with that of the quark, 

e2g2 

agg^^^i+i- (s) = o-o-^ (6) 

where (Jq = and the overall color factor = ^ arises from to the fact that only 

when the color of the quark matches the color of the antiquark can annihilation into a 
color-singlet, leptonic, final state take place. 

The quantity s introduced here, as well as t which will be introduced below are the 
partonic Mandelstam variables while s and t are these defined for the incoming hadrons. 
In general, the incoming quark and antiquark will have a spectrum of CM energies and it is 
more appropriate to consider the hadronic differential cross section In order to obtain 
this, one can start from the CM frame of the 2 hadrons. In this frame, the four momenta. 
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and P2 of the incoming partons may be written as 

/s 

Pi = ^ (a;i,0,0,xi) 

P2 = ^ (2;2,o,o, -X2) . 

The square of the parton CM energy s is related to the corresponding hadronic quantity 
by s = X1X2S. Folding in the parton distribution functions for the initial state quarks and 
antiquarks gives the hadronic differential cross section in terms of Xi and X2, 

= [^,ixuQ)T,{x2,Q) + {1 ^ 2}] (7) 

where is the parton density function of species q. The quantity Q is the factorization 
scale, usually taken to be the invariant mass, a/5. From beam symmetry, the substitution 
{1^2} is equivalent to simply multiplying the whole expression by 2. To obtain the 
hadronic differential cross section the transformation from Xi and X2 to s and y is necessary, 

S = SXiXo 

(8) 

where y is the rapidity of the pair. The transformation given in Eq ^ involves a Jacobian 
which reduces to the constant It is also necessary to integrate over all possible rapidity 
values, where this integration is usually done numerically. By knowing the value of s (14 
TeV for the LHC), the hadronic differential cross section can written as 

+yo 

da <5"o 2 /" 



ds 



^2 j ^^g2j^^j^^^^^^g)^y|j^_j^^^^^^g)^y|j^|^^2}] (9) 



where the boundaries ±?/o are determined from both Xi and X2 being constrained between 
and 1 so that ?/o = | In (f ) • It is sometimes useful to replace the transformation Jacobian 
with the equivalent expression, ^ = 

3.2.2 Partonic level: qq ^ J ^ l^l~ 

The differential cross-section function which PythiaS takes should describe the hard-process 
itself. It should not be given as the (hadronic) differential cross section and it should not 
include any parton density functions. This is because s is being determined separately 
for every generated event and since the evolution of the parton distributions is performed 
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internally by PythiaS. This is similar to the e~^e~ — ?■ 7 — )■ fi'^fi~ process considered 
previously where only the hard process was considered. Inserting the t dependency, the 
function given to PythiaS becomes 

di = l2"lFiVjT (l+4V.cos^) (10) 

where this is written in the CM frame of the incoming partons. The angle 6* is the polar 
angle in the CM frame between the incoming q and the outgoing /~ in contrast to the polar 
angle, 6, in the lab frame. 

3.2.3 A qq final state 

So far account has only been taken for the s-channel 7 exchange, since only the leptonic 
final state, different from the initial state, at only low energies has been considered. To 
generalize this Drell-Yan annihilation to any pair of fermions, gg — ?■ 7 — ?■ //, Eq |l^ must 
be modified. The cross section has to be multiplied by the appropriate color factor A'^^ 
and all the existing lepton indices must be substituted with corresponding fermion indices 
/ ^ /. The new overall color factor is A'^^ = 3 for / = g or A'^^ = 1 for / = /. Doing so, 
it is apparent that in the general di-fermion final state the contribution of the t-channel 
exchange should also be considered. This is since the di-fermion final state can consist 
of the same (annihilated) quark-antiquark pair, qq ff = (qq) sa_ync exchange 
can take place in either of the s-channel or the t-channel either with photon or gluon 
exchange as can be seen in Fig |^. Thus, a di-jet final state can also be observed. Indeed, 
the cross section should be much larger than for the di-lepton final state. However, the 
contribution of the s-channel 7 exchange to this di-jet final state is negligible since the 
electroweak interaction is much weaker than the strong interaction. Therefore, this channel 
will be dominated by the exchange of a gluon, either in the s- or the t-channel. The t- 
channel exchange is expected to be significantly larger than the corresponding s-channel 
exchange, in particular at forward angles. However, experimentally it depends on the 
transverse momentum cut, pt, applied on the outgoing quark, or jet, hadronic state. The 
t-channel exchange predominantly produces events with smaller pt- In addition, there are 
significantly more QCD diagrams involving gluons in the initial and final states which can 
lead to di-jet production. In practice, it is not possible to distinguish between the initial 
state and the final state sources. It is very difficult to distinguish between a quark and 
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Figure 2: The competing amplitudes for a production of a di-quark final state. The 
notations c, c', c, c' account for the correct color flow. The leading contributions come from 
the gluons exchange. 

gluon jet, unless for instance, it can be identified as a 6-jet by the presence of a displaced 
vertex. In principle it may be possible to suppress the t-channel contributions by a rejection 
based on the di-jet system px- Nevertheless, for simplicity, the remainder of this paper will 
address only di-lepton production from s-channel exchange. 



3.2.4 Including the Z° boson: qq j/Z° l^l at a/I > m^o 

With increasing energies, around and above the pole, the s-channel intermediate 
state and interference terms must be included. The helicity amplitude of the Z^ boson has 
to be added to the photon amplitude, egCi/s, to account for the correct interference effects. 



+ 



9\q9h 



^20 



F 



Here, the SM coupling constants |T^ of the Z° gauge boson to the involved fermions are. 



9x, 



sm cos Bw J 



li — er s\v? 9iY _ 

- — if A/ 



'1/2 



(12) 



^ sin 614/ cos Svi' '^f 

where the index / represents either the incoming quark or the outgoing lepton. Here the 
additional values for the Z^ boson mass, m^o, the weak isospin of the incoming quarks or 
the outgoing leptons, and the weak-mixing angle, 9w have been introduced. The sum 
z^^FF is the total Z^ decay width to all fermion-antifermion pairs denoted by FF, 

F 

where the partial decay width in Lowest Order (LO) is: 

N^aemmzo ^ 2 



6 



+ 



(13) 
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To account for radiative corrections, the relation = m^o sin^ 6w cos^ 6^ has been 

used, where is the muon decay constant and mw is the W boson mass. The corrected 
expression for the partial decay width is 

[Ilf - 2lleF sin^ Ow + I {ep sin^ Owf] • (14) 



37rv/2 

Finally, the differential cross section in its LO takes the form 

^^T^-l^-f 4t E_ Ej^^M.wr(i+«.v-n^. (15) 

In this context it is worth mentioning that self-formulated coupling constants were used 
in Eq |T^. These can be transformed into the PythiaS coupling constants, af,Vf,Lf,Rf 
given by the CoupEW class. The couplings used in MoSES can be modified without causing 
undesired behavior in PythiaS. Thus, they provide secure flexibility to probe non-SM 
phenomena that depend on them. The same considerations are also relevant for the widths. 

3.2.5 The angular decay asymmetry 

To expand the previous discussion it is useful to define two additional reference frames: 
(a) the colliding proton CM frame denoted by O (this frame is identical to the laboratory 
frame) and, (b) the rest frame of the di-lepton system denoted by O*. Neglecting higher 
order processes, the di-lepton system is, in general boosted along the beam axis. The 2;-axis 
is arbitrarily chosen as the direction of one of the beams, and it is then identical for O 
and O* frames. One of the primary observables for this process is the forward-backward 
asymmetry that can be extracted from the distribution of cos 6** in the O* frame. By 
definition, it is the cosine of the angle between the quark and the lepton directions in 
the O* frame. It should be noted that there is a sign ambiguity in the measurement of 
cos 6'*, since for a particular event, there is no information about whether the incoming 
quark comes from the positive or negative z directions. Instead, it is useful to consider the 
quantity cosO*^, where 9*^ is the angle between the di-lepton system boost (3 (relative to the 
O frame) and the lepton direction 

cos^^ = , Tr^-^ (16) 

^ \p*l\-\p\ 

where the boost vector is /3 = eIVe- - order to obtain p*i, the boost vector of the di- 
lepton system should be found and the transformation to the O* frame should be performed. 
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Neglecting higher order processes, the boost is confined to the \z\ direction, that is, /3 = /35 
and can be measured so that there is no sign ambiguity in determining cos 6'^. 

The next step is to calculate the forward-backward asymmetry Aj^. Here we define 

daf /ds — dab/ds daf /ds — da^/ds 



■^^^^^ daf /ds + dat/ds datot/ds 
where, in the absence of any detector cuts applied, the quantity daf/^/ds is given by 

1 



(17) 



dah 
ds 



da is) 



dsd cos 9 



-dcos 6* 



da 



f 



ds 



da (s) 



dsd cos 9* 



dcosO* . 



-1 

Since the integrand is given in terms of cos 6* it should be rewritten to enable re-classification 
of the forward and backward definitions. This can be done by separating into different ra- 
pidity contributions since y and /3 have the same sign. 



ds 



ds 



I dyj 



do-{s,y) 
dsd cos 9* 



+yo 



d cos 9* + f dy f 



d(T{s,y) 



dsd cos 9 



-dcosO* 



-yo 



(19) 



+yo 1 

-yo -1 



where yo is the rapidity kinematic limit (see Eq |^) and the form of Eq |T^ remains the same 
under the substitutions cr//b — cr^^^ and Aji, ^ A*^^. The integration over y in Eq ^ is 
performed asymmetrically and thus, the relation in Eq ^ cannot be used. Therefore, in the 
hadronic level the differential cross section is usually separated into two terms with respect 
to cos^*: Symmetric and Anti-symmetric denoted by S and A. This enables the separate 
contributions that form the forward-backward asymmetry to be identified, 

dyds7,^,e* ^ 5Z ^) ^« (^) (1 + + (y^ ^) ^1 (^) 2 cos r ] (20) 

where the Symmetric and Anti-symmetric combinations Gq^^~^^ that involve the parton 



density functions of the colliding hadrons are defined 



X1X2 



J^q ( Xi, a/I) ( X2, V^) ± J^q \ Xi, J^q ( X2, 



(21) 



where xi and X2 can be written in terms of y and s so Gq is in fact a function of y and s. 
The terms Sq and Aq can be realized from the helicity amplitude: 



EE 

Xq A; 

^ E 

Xq—Xl 



MXqX, 
XqXl 



E 



M 



XqXl 



(22) 
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where Ma^A; from Eq ^ is sensitive to the couphngs of the Z° boson to the different famihes 

of SM fermions. Using Equations ^ and it can be seen that under these definitions 

the forward backward-asymmetry expression reduces to 

+yo 
E A, (s) f dyGf {y, 5) - E (s) f dyGf {y, s) 

A% is) ^ ^ ^ ^ . (23) 

E S, (s) J dyGS {y, s) 
1 -yo 

In cases where the integration over cos 6* is performed in a smaller effective interval {[—k, k] 
where < k < 1) due to some kinematic cuts, then the substitutions j dy ^ J dyk"^ and 
j dy ^ J j dy {k + ^ j can be made in the numerator and the denominator respectively. 
These substitutions will be relevant in the following discussion where the ATLAS detector 
cuts are introduced. 

Thus, the origin of the asymmetry is understood and can be summarized by simply 
considering the cos 9* distribution 

^ \[l + cos^ei]+A%co,ei. (24) 



NdcosOl 8 

The forward-backward asymmetry coefficient depends on the couplings of the fermions 
to the boson and is thus, sensitive to sin^ 9w and it can be extracted by fitting Eq. ^ 
to the data. However, in the case of high statistics and where there are no kinematic cuts, 
a direct measurement could be done by simply counting the forward and backward events 
Nj^f^. These are manifestly given relative to the boost direction if cos 6'^ was first used 

^ _Nl-N[ 

~ l^fT^- ^^^^ 

A detailed discussion of the polar angle distribution and the asymmetry is given in the 
following sections including the effects of Initial State Radiation (ISR). 

3.2.6 Validation results with j/Z^ 

Using the expression given in Eq |l^ events were generated simulating colliding proton 
beams at the CM energy of 14 TeV with no cuts applied, and all the parton-level switches 
were turned off. Considering only the di-muon final state, all other decay modes of the 
'y/Z^ were also turned off. The same validation procedure described previously is repeated, 
combined with the comparison of the hadronic cross section shapes. As a second validation 
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stage, the code PythiaS uses internally was copied and was plugged-in as if it was an 
external process. In the following this is referred to as the semi-external validation. 

The differential cross section shapes and the distributions of the cos 6*. are shown in 
Figures ^a,c and ^a,c for the external and semi-external process implementations respec- 
tively. Corresponding bin-by-bin comparisons are shown in Figures ^b,d and ^b,d. The 
numerical results are summarized in Table ^ Note that since in this section, the gener- 
ation is stopped at the hard process level, the asymmetry results are likely to be slightly 
modified by initial state radiation and will be discussed in a later section, but this should 
not, in practice, modify the agreement between the results of the three implementations. 
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Figure 3: The results obtained with MoSES external code (sohd) and PythiaS inter- 
nal code (dashed). The \/I distributions and their comparison in (a) and (b) and, the 
normalized cos 6*^ distributions around the peak and their comparison in (c) and (d) . 

The agreement between the total cross section values in Table ^ is to within ~3-sigma, 
where most of this difference might be due to the phase space sampling. The forward- 
backward asymmetry values in Table 0, extracted from the three cos 6*^ distributions, agree 
to within ~0.1-sigma. Figures ^b,d and ^b,d illustrate the agreement of differential dis- 
tributions to within the respective statistical uncertainties. A comparison between the 

2 

internal and external histograms yields j^^p = 1.04 for the di-lepton mass distribution and 
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Figure 4: The results obtained with MoSES semi- external code (sohd) and PythiaS 
internal code (dashed). The distributions description is identical to Fig |. 

Table 2: Numerical comparisons. The cross section statistics as given by PythiaS and 
the A^j^ value around the resonance. This is the output for the three IM samples of 
pp ^ ^ / fi~^fi~ events with 14 TeV CM energy. The lower mass cut-off is at 20 GeV. 



Source code Estimated crtotai[nb] 



Events in 



^ ptot 



Internal 
Semi- External 
External 



2.624 ±0.001 

2.618 ±0.001 

2.619 ±0.001 



0.0445 ±0.0018 
0.0459 ±0.0018 
0.0469 ±0.0018 



308160 
309381 
309186 



2 

= 1.03 for the cos 6'^ distribution. Likewise, the comparison test between the inter- 

2 

nal and semi-external histograms yielded = 0.94 for the di-lepton mass distribution 



2 

and = 1.12 for the cos^|, distribution. 
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4 The new KK boson exchange process 

The Drell-Yan KK process, pp — )■ '-f*/Z* l^l^X, is similar to those discussed in the 
previous section and has not yet been implemented as an intrinsic component of any public 
Monte Carlo generator. After following the validation procedure for the formalism above, 
it can be used to introduce the new specific KK process within the theoretical framework 
given below. Also given below is a brief discussion on other possible heavy gauge bosons 
that can be produced in the LHC and share the same final di-lepton state. 



4.1 An overview of the theoretical framework 
4.1.1 The Kaluza-KIein model 

The observable world consists of at least the three known spatial dimensions, and one of 
time, but it is possible that there are additional spatial dimensions that are not directly 
observable. When considering these potential extra dimensions, it is useful to distinguish 
between parallel and transverse dimensions with respect to the 3d world. The size of all 
parallel dimensions should be constrained to be no larger than ~ TeV~^ (10~^^ m) in 
order to be unobservable at present energies. The size of transverse dimensions remain 
unrestricted with much weaker experimental bounds. One popular possibility is 



one extra parallel dimension, compactified on the /Z2 orbifold where gauge fields have 
KK excitations in this extra parallel dimension but where fermions are localized on the "id 
brane and have no KK excitations. KK states of gauge bosons are then singly produced 
as new resonances that may be observed experimentally. This description leads to equally 
spaced KK states of gauge fields with masses given by 

ml = ml + ^ (26) 

where m„ is the higher dimensional mass and R is the radius of compactification. The mode 
= is identified with the 4d (SM) state, while the higher modes have the same quantum 
numbers as the lowest, but with increasing mass. Within this model, the couplings to 
fermions of these excited modes are larger than the known couplings of the zero modes 
(SM bosons) by a factor of a/2 due to the normalization of the KK excitations []7|,P^ 



In this paper, only in the excitations of the Z° and the photon shall be considered. If 
such KK states exist, they have to exceed the lower bounds on their mass, based on their 
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indirect effects associated with their tower exchange. These bounds rely upon a number 
of additional assumptions, notably, that the effect of KK exchange is the only new physics 
beyond the SM. In the 5d case, a global fit to the precision electroweak data including the 
contributions from KK gauge interactions yields >4 TeV ||T3|, |20|-E3] . The reader is 



referred to the appendix for a derivation of the KK tower for a 5d real massless scalar field. 
In the following discussion, the notation m* = is used. 

4.1.2 Additional heavy Z-like boson production 

A similar deviation from the SM cross section coming from KK excitations of 7/2'°, may 
be observed also in models with an additional heavy gauge boson |]|,^,[12|. Several Grand 
Unified Theories (GUTs) postulate that the SU{3), SU{2) and U{1) symmetry groups 
of the SM have a common origin as sub-groups of some larger symmetry group G. It is 
supposed that at large energy scales, this symmetry is valid but below some critical energy 
scale, G, it is spontaneously broken. This kind of GUT predicts at least one additional 
gauge boson after the symmetry is broken to the SM. In these models, the extra neutral 
gauge bosons are usually denoted by Z'. 

For comparison with the KK model, one specific model that is featuring at least one Z' 
boson (light enough to be detected at the LHC) was considered. These models can come 
from the breaking of the Eq group which is a popular candidate to GUT symmetry. The 
Z' can be observed as a peak in the di-lepton mass distribution above a small background. 



and the LHC discovery potential for that is reasonably high and well known p3[. If a 
resonance is observed at the Z' or 7*/Z* hypothetical mass at the LHC, a discrimination 
mechanism between these two candidates would be required. In specific circumstances. 



this discrimination is possible |]T2|,^, however, for both the generation of the KK events 



and the discrimination mechanism, a somewhat different approach is considered in this 
paper. Out of the possible ways to break the Eq group, the following shall be considered 
(z) Z'^:Ee^ SO{10) x f/(l)^, (zz), Z'^:E^^ 50(10) x U{1)^, ^ SU{h) x U{1)^ x f/(l)^ 
and (Hi) Z'^ : Eq ^ SM x f/(l).;y. The different couplings of the new gauge field to the 
SM fermions within these models can be found elsewhere [0. Another scenario often 



introduced is the Zgj^ where the new boson has the same couplings as the Z^ but with 
different mass and width. One should note that there is no theoretical justification for the 
choice of SM-like couplings for the Z'. However, this is the more experimentally challenging 
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case since if there will be a resonance in the KK mass, then the case where the couplings of 
the Z' are SM-like is practically the most difficult to distinguish from the KK resonance. 
Therefore, this will be the choice for the following discussion, where this specific model 
shall be denoted by Z'^y^. In this paper, the mass of these new bosons are taken to be 
4 TeV, same as the mass of the KK bosons. 



4.2 The KK implementation in the LHC scenario 

As mentioned in the introduction, the LHC represents a new frontier in the search for heavy 
resonances. Each of the resonances in the KK tower discussed here could be produced by 
a similar mechanism to the light SM bosons. From an experimental stand point, most of 
the considerations described in the previous sections hold also for the KK case. 

In the parton level, the process gg — > 7*/^* — Z"*"/^ can be expressed in terms of the 
following differential cross section 



rfa(g,COSr) ^ 2 OL^S^ 



(n) 
n=0 



l + 4:Xg\icos9*f (27) 



where the complete amplitude consists of the SM term, exactly the term given in Eq [IT], 
plus an infinite KK tower of excitations with increasing mass. 



00 00 



n=0 n=l 

where Mi'^l = M\ \, and where each contribution for n > 1 can be written as 

(n) (n) 

is) ^ + ^ . (29) 



Recalling Eq |2B|, the n*'^ KK excitation masses m!'^} and rrS^} are given by 



m^i = y/mlo + (n ■ m*) 



*^2 



in) if 

m^, = n ■ m 



(30) 



As discussed previously, the current limits on the m* value are approaching 4 TeV and the 
LHC is expected to enable the expansion of the the search region. Practically, the mass 
m* =4 TeV is the value taken for the KK amplitudes as well as arbitrarily choosing an 
upper limit of n = 100 which is large enough so that in practice, higher excitations do not 
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play a significant role within the accessible LHC energy range. Even though all the KK 
excitations higher than the first are beyond the reach of the LHC, their presence still affects 
the accessible LHC energy range due to interference contributions which are non-negligible 
even at energies far below their masses. As mentioned, the couplings of the excited KK 
states to fermions are larger than the SM ones, Eq |12|, by a factor of ^2; 



(n) 



gx^ if n = 
\/2 ■ gx^ otherwise 

(31) 

6/ if n = 
\/2 ■ e f otherwise 

where as in the previous section, the index / can either represent the incoming quark or 
the outgoing lepton. The sums X^^^i/ *^pp are the total 7* and Z* decay widths to all 



(n) 



fermion-antifermion pairs denoted by FF. Considering the SM terms from Equations 113 



and [T^, these forms are affected by the a/2 factor introduced in Eq In addition, with 
respect to the SM terms, a single power of the mass appearing in Eq |1^, is replaced with 
the nominal mass (See Eq RH), starting from Eq for the Z* and from Eq 113] for the 7*, 



i Z*^FF - ZO^FF ■ \ ^^(n) 



1 if n = 



'r*^FF 



2—^ otherwise 

(32) 

NPa,mm["J I if n = 

4e|n otherwise 



This process can be realized as shown in Fig |^. 

Form the MC event generator stand point, in order to account for the entire KK tower 
- namely a large enough n - the amplitude itself should be coded as a complex expression 
so that all the interference terms will emerge. For instance, in the 'j/Z^ case in PythiaS, 
the amplitudes are pre-calculated and separated into the three different contributions (7, 
Z° and interference). Hence it is clear that in case of large number of diagrams, there will 
be many consequent interference terms and it is difficult to calculate and code. 
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[n = 0) {n = 1) 

n (0) (1) * (1) / 2 ; 12 

nu^ — m^o = m^o m:^: = m rrv^i = y m^o + rn*^ 

Q /+ 9 1+ Q 1+ Q /+ 

Q 1-1 1-1 1-1 l- 

Figure 5: The KK excitations tower of the gauge bosons '~f/Z^ starting from the O*'' SM 
state. Note the couphngs and masses (which affect also the widths) of each level. 

4.3 Analytic results 

In the following discussion, the general behavior of the cross section is studied under the 
LHC conditions, within the kinematic range of interest defined by several selection criteria. 

Specifically, the di-muon final state, pp — )■ '-/* /Z* — )■ fi~^fi-, is chosen for this preliminary 
study. The range of interest is determined by > 10 GeV and |?7| < 2.5 which corresponds 
approximately to the ATLAS trigger and detector acceptance. To complete this definition, 
another threshold on the di-muon invariant mass is considered: ^/I > 1 TeV. 

The restricted range implies that the hadronic distributions given in Equations ^ ^ ^ 
and ^ should be modified since they were obtained by integrating over cos^* within 
I cos 6** I < 1, whereas the effect of the cuts is restricting it to | cos 6^*1 < |cos^i!Qax| given by 



^cut \ 2 



IcosCaxI = min { |tanh (r/™* - IyI) h \/ 1 " ( ^ ) } (33) 



where r] = 2.5, p'^ = 10 GeV, and y is the rapidity of the lepton pair. Due to these cuts, 
y and \/§ are consequently confined to \y\ < 2.5 and > 2p^^. Using the same analytic 
approach, the two expressions in Eq ^ can be inverted to obtain the effective cuts on s. 

The y/l distributions around the first excitation are shown in Fig ^ where a strong 
destructive interference is clearly seen between 1 and 2 TeV for the KK line. The distri- 
butions of the SM and the Zgj^ model discussed in the previous section are also shown for 
comparison. As mentioned in the introduction, the shapes from Fig ^ could enable the 
discrimination between the models already at C = 100 fb~^. 

An estimation for the number of events that are expected to be measured by the ATLAS 
detector in the range of interest is necessary for generating the MC pseudo-data samples 
with a realistic size. The mass range considered is, 1 < < 6 TeV with integrated 
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luminosity of C =100 fb^^. After applying the selection cuts a total number of ~400 KK 
events are expected in that overall mass range at the respective luminosity and ~190 KK 
events in the KK peak area, 2 < < 5 TeV, see Table ^. The number of events in the 
overall mass range of interest, 1 < < 6 TeV, will be used as an input for the simulation. 

The analytical function for the forward-backward asymmetry relative to the boost di- 
rection, denoted by A^^, is shown in Fig|^. Large differences between the forward-backward 
asymmetry for the Zgj^ and the KK models around the KK resonance are apparent. De- 
pending on the integrated luminosity, these expected differences, together with the obser- 
vation of a significant peak, might enable the discrimination between the three models. As 
illustrated in Table ^ the statistics in this range of masses is expected to be poor and as 
such, calculating the asymmetry in a wider region might be necessary. To give a prediction 
for y4yj, in 2 < < 5 TeV, the theoretical function is averaged over y/I, 

5^/|4WSx(/|.f'. (34) 

The predicted values for this averaged asymmetry, A^^, in LO are summarized in Table ^ 

Table 3: The expected number of events in C=100 fb^^ for tie tliree models, SM, Z'^-^ and 
KK for 1 < y/I < 6 TeV and for 2 < < 5 TeV within the described kinematic region. 



Model 


Events in 1 < < 6 TeV 


Events in 2 < < 5 TeV 


SM 


~ 480 


~ 15 


7' 


~ 460 


~ 30 


KK 


~ 400 


~ 190 





Table 4: The expected values of in LO for the three models SM, Z'^^ and KK around 
the KK resonance 2 < \/3 < 5 TeV within the described kinematic region. 



Model 


A'j^m2<VI<5 TeV 


SM 


0.325 


7' 


0.090 


KK 


0.308 
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Figure 6: The invariant-mass distributions of the three models, KK (sohd), Zgj^ (dotted) 
and SM (dash-dot) within the kinematic region described in the text. 
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Figure 7: The forward- backward asymmetry within the kinematic region described in the 
text, corresponding to the three models, KK (solid), (dotted) and SM (dash-dot). 
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4.4 Strategy of the experimental kinematics study 

This sub-section presents some preliminary results obtained using MoSES for the KK 
process with the di-muon final state, pp — ?■ '-/* /Z* — )■ + X, at the generator level. 

The results from the KK process are presented and compared with the SM and the Zgjyj 
results (with the same final state) with an emphasis on the discrimination between the 
possible KK or Z'^y^ signals. 

The realization of the processes within the LHC and ATLAS kinematic regions enables 
observables to be studied in a realistic regime at the generator level. The next stage is to 
embed the new subprocesses within fully simulated events including the effects of initial 
and final state radiation (ISR and FSR), parton showering, hadronization, proton remnant 
fragmentation, particle decay etc. All these effects are available within PythiaS. Because 
of factorization, only the hard subprocess needs to be generated, and integrated within 
PythiaS. Therefore, PythiaS itself is responsible for the subsequent forward evolution of 
the event from the time and energy scale of the provided hard-process. It is also responsible 
for the backwards evolution to the initial conditions of the interacting protons. 

The analysis starts from the measured di-lepton final state, namely, the 4-momenta of 
the charged leptons p'( and pj given in the colliding protons CM frame, denoted by O. The 
sum of these 4-momenta measured in the O frame is 

g^=pr+pr (35) 

It is useful to define the di-lepton squared invariant mass, = QyQ'^ = s and its 
transverse momentum, Qt = [pj +p},p'i +Pj)- The complementary components, = 

{Q'^ ± Q^), can also be defined It is useful to write the rapidity of the di-lepton system 
l/Q = I In (^Q^^ ■ For the following discussion the masses of the leptons and partons are 
small with respect to Q = \f§ and are thus neglected, as are the proton masses which are 
small with respect to ^/s. 

In the LO Drell-Yan picture, both original partons have momentum only along the 
original z direction of the O frame. Correspondingly, the di-lepton system CM can have 
momentum only along the same direction. When considering higher orders in the Drell-Yan 
picture, if one of the original partons emits an ISR gluon it can obtain momentum also in 
the transverse direction. The intermediate state and the di-lepton system will then have a 
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momentum component in the transverse direction. In general this will modify the angular 
distribution of the outgoing leptons with respect to that discussed earlier. In addition, 
it is also possible for one of the outgoing leptons to radiate an FSR photon. If it is not 
detected and corrected for in the lepton kinematics, then the resulting invariant mass of 
the di-lepton system will be shifted towards lower values. The effects of the ISR and the 
FSR can be demonstrated in Fig ||. 




"^sdi-ii^HGeV] 

Figure 8: The Qt vs. \f§ distribution of the exchanged KK bosons. The separate contri- 
butions of ISR and FSR are clearly seen (See text). 

The unique 2d distribution seen in Fig ^ is the first comprehensive result of the entire 
Moses and Pythia8 framework that is being presented. It illustrates several aspects of the 
KK process which can be seen in it very clearly. To obtain this plot, a sample of 10^ KK 
full events was generated using MoSES and Pythia8, including all steps of the complete 
event generation. The events were selected such that they are within the kinematic range 
of interest introduced in the previous sub-section. The KK resonance can be seen around 
~4 TeV with two distinct Qx tails, one vertical which is associated with the ISR 
effect and the other one diagonal which is associated with the FSR effect. This picture, of 
initial and final state radiation, piles up more complication as can be clarified in the next 
discussion. 



25 



Thus, from the Qt distribution of the di-lepton system (with respect to O frame), it 
can be clearly seen that the di-lepton CM has a non-negligible momentum in the transverse 
direction due to the radiation of the initial state gluon or final state photon. The tails seen 
in the Qt distribution in Fig |^ can be rather large, causing once more a problem for the 
way cos 6** is defined. 

Following the results seen in Fig ^, if the FSR switch in PythiaS is turned off then the 
related Qt diagonal tail vanishes, as expected, and no corresponding shift in the invariant 
mass occurs. In addition, turning off the ISR as well, then the vertical Qt tail does 
not extend significantly above at Qt ~10 GeV rather than the ~1 TeV shown in Fig ^ 
However, further consideration of these phenomena suggests that they may not have a 
significant effect on the overall kinematic distributions. 

For di-muon events where the FSR photon has sufficiently high transverse energy, the 
measured muon might fail the pt > 10 GeV cut and this would affect the distribution. 
However, in this case, the photon itself should be sufficiently well separated from the muon 
so that it can be detected independently. On the other hand, for di-electron events with a 
photon being sufficiently collinear with the electron, the calorimeter cluster would include 
both the electron and photon and the reconstructed energy would correspond well to the 
initial electron energy before radiation. In any real study, the effect of final state photon 
radiation can be simulated and the effect on the trigger and selection efficiency can be 
investigated. However, this case will not be discussed further here. The case of ISR gluon 
requires more attention and will be discussed in the next section. 

4.4.1 Handling with the ISR effect 

To consistently handle the effects of ISR, since the boost direction no longer coincides with 
the beam axis, it is not clear which axis should be used for the measurement of cos^*. 
There are two approaches to deal with that problem; 

• minimize the ISR and FSR effects by applying a cut on large Qt's {ie, Qt < Qt^^ ~ 
(9(100 GeV)), where for most of these events, the di-muon system momentum will be 
along the beam axis. In this case, the LO formalism where there is no ISR effect is 
valid and can be used. However this will reduce the number of events in the sample. 

• analyze the events while taking into account the ISR, by rotating the di-lepton rest 
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frame by some angle so that the 2-axis will be slightly modified. A common choice 
is the Collins-Soper (CS) reference frame [^5| -p7|, denoted here by O'. It can be 



shown that this choice, minimizes the contribution from longitudinally polarized 
'~f*/Z* which can be produced due to the gluon ISR and thus, can affect the an- 
gular distribution. 

Note that in the O' frame, cos 6' again has the sign ambiguity previously discussed. This 
ambiguity results from our arbitrary selection of the original z direction in the O frame. 
This implies that a reclassification of cos 9' can be performed with respect to the rapidity 
sign such that if yq < then cos 6'' — > — cos^'. This is possible since in most of the 
events, the CS rotation angle will be small enough so the choice of the z' axis will point 
approximately in the direction of one of the protons (p'p^ or p'^^)- 

The lepton cos^' distribution to first order in 0{as) must have the form p6 



dN 



NdcosO' 8 



1 + ^Aq + ^4 cos 6'' + M - ^Aq ) cos^ 9' 



(36) 



where the coefficients Aq and A4 can be functions of the kinematic variables s, s, yq 
and Qti and where for each event cos^' should be calculated relative to the sign of yq. 
Finally, the forward-backward asymmetry in the CS O' frame can be realized from the last 
expression as A'^f = ^A^ and, as expected, now depends on both s and Qt- The resulting 
forward-backward asymmetry should be unaffected by the gluon ISR and the consequent 
transformation to the CS frame, so the only affected part is the symmetric one that involves 
Aq. This coefficient is relevant for the examination of the complete cos^' distribution. 

The angular distribution given in Eq ^ corresponds to an exchanged particle with 
spin-1 and in that context, the cos^' distribution of a spin-2 particle such as the graviton. 



exhibits a completely different behavior |114|, p!5|, |28| . As mentioned in the introduction, a 
spin classification should be also done. 
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4.4.2 The motivation for extracting the angular coefficients 

The motivation for extracting the A4 coefficient is the insight it can provide on the couplings 
of the new heavy gauge bosons to the SM fermion fields. Apart from that, there is a 
possibility that this coefficient may support the discrimination between the two spin-1 
models, as well as a spin-2 model. However, it is clear that extracting the coefficients 
from the pseudo-data and comparing with the MC reference coefficients, can give no more 
discrimination than comparing the distributions themselves. 

As illustrated in Table § and in Fig ^ the two non-SM forward-backward asymmetry 
coefficients differ significantly for any mass in the range 2 < < 5 TeV. However, the 
values given in Table considered as an analytic results, are averaged and also not strictly 
accurate for the following analysis since the ISR (and FSR) effect was not included. In 
the case of ISR, the measured asymmetry is expected to slightly change and in addition, 
there is another coefficient, the symmetric Aq (vanishing in the case of no ISR). Thus, 
the best estimation for these coefficients, under the influence of ISR, can come from a 
high statistics MC reference sample^ (for each model). The coefficients can be extracted 
by fitting the cos 6*' distribution (given in Eg [36|) corresponding to the high statistics MC 
reference sample and this is done also for the realistic statistics (e.g. pseudo-data with 
C = 100 fb~^) to quantify how well it might work. In that sense, it should be stressed 
that although the Aq coefficient is relevant for the complete cos 6'' distribution study, its 
role here is no more than to technically improve the fit whose primary objective is the 
forward-backward asymmetry coefficient. The cuts on px, rj and y/S, introduced in the 
previous section, are applied on all the samples. 

4.4.3 The discrimination between the spin-1 models 

As mentioned in the introduction section, having the large MC reference samples and 
the £=100(500) fb~^ pseudo-data samples, the comparison can be performed using the 
Kolmogorov test for the \/S and the cos 6'' distributions. Each distribution exhibits a 
special characteristic behavior for the examined model. In that context, the Kolmogorov 
test can also provide good sensitivity for the determination of the spin of the exchanged 
particle. 

^high statistics compared to the expected between 1 TeV and 6 TeV, see Table ^ 
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4.4.4 The strategy 

The strategy used to classify the observed signal, can be summarized in four points: 

1. searching for a significant resonance in the di-muon distribution above the small 
expected SM background. 

2. comparing the cos 6'' distribution around the resonance to a spin-1 and spin-2 reso- 
nance distributions. This step is not shown in this paper]]. 

assuming a spin-1 resonance, 

3. comparing the y/I and cos 6'' distributions to all three MC reference distributions. 

4. fitting the cos 9' distribution from Eq ^ to the data. 

In the next sub-section, a collection of results obtained from the preliminary analysis at 
the generator level are gathered and shown in detail. 

4.5 A collection of the results 

In the first part of this subsection, a detailed discussion is given on the Maximum Likelihood 
(ML) fit results for the cos^' distributions. In the second part, the Kolmogorov test 
results for the cos 6' and y/I distributions are summarized with detailed comments on the 
possibilities of discriminations. 

First, let the pt, rj and cos 6'' distributions, for the three models, be gathered and 
shown at high and low statistics. The kinematic distributions of the large MC reference 
samples are shown in Fig ^ where the distributions of the pseudo-data samples can be seen 
in Fig ^ and |Tl| respectively. 

In Fig ^ the sizes (luminosities) of the SM and the Zgp^ MC reference samples are 
normalized to the arbitrarily chosen size (luminosity) of the KK reference sample (10^ KK 
events). In Figures |10] and [11] the sizes of all three samples are determined by the two LHC 
luminosities, ie £=500 fb~^ and £=100 fb"^. These three figures are arranged in the same 
format so their properties are identical. 



Optionally, the azimuthal angle distributions can also be used for this purpose 



29 



woo 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000 



(a) 



10 













-■ SM - 






KK 1 




1.. 


L 7' - 


















, \^f^n, -; 









500 1000 1500 2000 2500 3000 3500 4000 



(b) 



-10 12 3 



(c) 




-O.S -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 
cosQ ' (Collins— Soper frame) 



(d) 



Figure 9: Kinematic distributions of the MC reference samples for the KK (sohd), the SM 
(dash-dot) and the Zgjyj (dotted) models, (a) The di-muon a/S distribution, (b) the muons 
Pt distribution, (c) the muons normalized rj distribution and, (d) the normalized muons 
cos 6' distribution. 
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Figure 10: Kinematic distributions of the C=500 fb ^ samples. See Figl 
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Figure 11: Kinematic distributions of the C=100 fb ^ samples. See Figl 
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4.5.1 The ML fit for the spin-1 cos 6' distribution 

The ML fit procedure, on one hand, should be unbiased by the influence of the speciflc 
kinematic range of interest and by the choice of a speciflc model. On the other hand, due 
for poor statistics one will want to use all events that are in the kinematic region. In other 
words, one is left with two possibilities: 

1. performing the flt in the entire allowed range of cos 6'' in order not to lose events at 
the chopped edges of the distribution^. The described kinematic region along with 
the choice of a certain model out of the three, will give a biased flt result. This bias 
will appear due to the fact that the angular acceptance, (which should be explicitly 
modeled in the flt), depends on the A^^ and thus, it depends on the choice of the 
model. Therefore, a detailed study of the dependence of the flt result for A^^ on its 
expected value should be performed. This study should be done for various values of 

so eventually one can understand this dependence and correct the flt result. 

2. performing the flt in a smaller range of cos 6'' where there are no edge effects due to 
the limited kinematic range of interest. This way part of the event sample is lost but 
there is no need to introduce a prior correction to compensate the bias tendency. 

Using the event-by-event ML method, the second option is chosen to flt the CS distribution 
(from Eq [36| ) to the large MC reference sample and the two pseudo-data samples. The 
range of the flt is limited to |cos^'| < 0.85 and that enables to perform the flt without any 
bias. The ML flt results for the Aq and A^^ coefficients of the MC reference samples 
and the pseudo-data (£=100(500) fb~^) are summarized in Table ^ The fitted curves of 



the MC reference samples are shown in Figures 0, and |14 



see the chopped edges of the cos 9' distribution in Fig ^ |3| and |lj 
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Table 5: A summary of the ML fit results for the cos 6' distributions in the CS O' frame 
around the KK resonance 2 < < 5 TeV and within the interval Icos^'l < 0.85. 







MC 


£=500 fb^i 


£=100 fb-i 




Model 


reference 


pseudo-data 


pseudo-data 




SM 


0.3476 ± 0.0053 


0.37 ±0.11 


0.48 ±0.25 




7' 

SM 


0.0900 ± 0.0042 


0.117 ±0.083 


0.24 ±0.18 




KK 


0.2958 ±0.0015 


0.216 ±0.035 


0.304 ±0.078 




SM 


0.00 ±0.02 


0.09 ±0.39 


0.62 ± 1.40 


Ao 


7' 


-0.025 ±0.015 


-0.53 ±0.26 


1.06 ±0.86 




KK 


0.0075 ±0.0056 


0.00 ±0.13 


0.00 ±0.28 




SM 


30145 


73 


13 


Events 


7' 

^SM 


57053 


142 


31 




KK 


413391 


798 


150 




cos^ ' ( Collins -Soper frame) 

Figure 12: The SM cos 9' distribution in the range 2 < < 5 TeV vs. the ML fit. 
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Figure 13: The Z^y^ cos 9' distribution in the range 2 < \/§ < 5 TeV vs. the ML ht. 
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Figure 14: The KK cos 6' distribution in the range 2 < < 5 TeV vs. the ML fit. 
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4.5.2 The Kolmogorov test for the cos 6' and \/l distributions 

The returned value of the Kolmogorov test is its probability, ie, a value much less than one 
means NOT compatible 0. The unbinned Kolmogorov test is chosen for two reasons; first, 
at low statistics it is usually better to perform an unbinned analysis and second, it does 
not add more arbitrary systematics that have to be studied due to sensitivity to arbitrary 
binning choices. The unbinned test results are summarized in Tables ^ and |^. 



Table 6: A summary of the unbinned Kolmogorov test results for the cos 6' distributions in 
the CS O' frame within the interval 2 < ^/I < 5 TeV and in the full cos 6' interval ([-1,1]). 



MC ref 


Luminosity 


T^Events of 


SM 


7' 


KK 


model 


of pseudo-data 


pseudo-data 


pseudo-data 


pseudo-data 


pseudo-data 




£=500 fb-i 


82 


0.97 


0.002 


0.008 


SM 


£=100 fb^i 












16 


0.887 


0.66 


0.89 




£=500 fb^i 


160 


0.018 


0.811 







£=100 fb-i 










35 


0.458 


0.458 


0.013 




£=500 fb'i 


971 


0.724 


0.011 


0.18 


KK 














£=100 fb-i 


181 


0.933 


0.822 


0.999 



Table 7: The unbinned Kolmogorov test for the y/I distributions within 1 < \f§ < 6 TeV. 



MC ref Luminosity ^^Events of 
model of pseudo-data pseudo-data 


SM Z'sM KK 
pseudo-data pseudo-data pseudo-data 


£=500 fb^i 2400 

SM 

£=100 fb-i 480 


0.254 0.002 
0.887 0.006 


£=500 fb^i 2300 

7' 

£=100 fb-i 460 


0.017 0.663 
0.28 0.016 


£=500 fb-i 1980 

KK 

£=100 fb^i 400 


0.308 
0.108 



In light of the results seen in Tables ^ and |^, it should be mentioned that a binned 
Kolmogorov test can provide better results {ie, the models will be more self-compatible 
and less inter-compatible), depending on the arbitrary choice of the binning. 
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4.5.3 A summary of the overall procedure 

1. Looking at Table ^ and comparing the values of Aq and Aft around the expected 
resonance, it can be concluded that: 

• Based on the fit results for the Aq^ coefficients, for the three models (in the CS 
O' frame), they turn out to be consistent with zero. 

• At £=100 fb^^, the values of A^^ for the three models are compatible with their 
MC reference estimations within less than one sigma. 

• At £=500 fb~^, the values of A^^ for the SM and the Zgp^ models are compatible 
with their MC reference estimations within less than one sigma, where the KK 
model is compatible only within 2.3 sigma. 

• The errors for A'^^ in both luminosities are too high. Hence, the sensitivity 
for probing the couplings using the measured A^^ coefficients requires larger 
samples. 

2. From the unbinned Kolmogorov tests in Table |]for the cos 6'' distributions, it can be 
concluded that: 

• There is clear distinction between the KK model (pseudo-data) and the Zgp^ 
model (MC ref) at £=500 fb^^ as well as at £=100 fb^^ 

• The Zgj^ model (pseudo-data) and the KK model (MC ref) are compatible at 
£=100 fb^^. However, the KK pseudo-data represents better its MC ref. 

3. From the unbinned Kolmogorov tests in Table |^ for the -y/S distributions, it can be 
concluded that: 

• As expected, there is a clear compatibility between all three models (pseudo- 
data) to the data simulated with the equivalent at £=500 fb~^ and even at 
£=100 fb^^ except for the Zg^^ model at £=100 fb~^ which is very small 
(-0.016) 

• There is a significant distinction between the KK model and the Zgj^ model at 
both luminosity values. 
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The important conclusion is that assuming an observed spin-1 resonance at \/§ ~ 4 TeV, 
above the SM background, it will be possible to distinguish between the KK and Z'^y^ 
models based on the Kolmogorov test for the a/I distributions, already at £=100 fb~^. 
For the cos 6' distributions, the Kolmogorov test and the ML fit for the forward-backward 
asymmetry measurement can be very important in providing supportive information for 
such a discrimination. This conclusion is valid already at £=100 fb~^ but it is much 
stronger for £=500 fb^^. In addition, it is clear that for both luminosities, the measurement 
of is too coarse for placing a precision statement about the new exotic couplings and 
it is only at higher integrated luminosities where a more sensitive study can be performed. 



5 Conclusions and Outlook 

In this paper the scenario of an observed resonance around 4 TeV, arising from the mea- 
surement of di-muon events is discussed. It is shown that for the LHC, it will be possible 
to discriminate between the specific KK model and the Z'^y^ model described in this paper 
assuming collisions in the design energy ^/s = 14 TeV and assuming integrated luminos- 
ity of £=100 fb~^. This statement relies both on the measurement of the \/§ and cos 6'' 
distributions and the measurement of the forward-backward asymmetry, applying as few 
cuts as necessary to deal with ISR. One should keep in mind that this conclusion can 
not be complete without passing the generated events through the full ATLAS detector 
simulation. In that context, it should also be commented that: 

• It is mandatory to apply a full detector simulation in order to treat the Kolmogorov 
unbinned (single) test as a single experiment. Furthermore, one needs to repeat 
these experiments with large number of generated pseudo-data samples to quantify 
the sensitivity of the test. 

• The comparison should also be done with a spin-2 model (namely, the RS graviton). 
In that context, the azimuthal angle distribution can be used as well [|1^]. 



A measurement of both di-muon and di-electron events will double the statistics 
shown here. 

Since the number of events is linear to the integrated luminosity, then it is clear that 
a sample much smaller than £=100 fb~^ will not be sufficient for any discrimination. 
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• In the first year of the LHC operation, it is expected that the colhsions will take 
place at y/s =7-10 TeV CM energy and the luminosity is not expected to exceed 0.2- 
0.3 fb~^. In this case it will not be possible to discover in the first year a resonance 
at 4 TeV related to any of the models discussed in this paper. 

• For the possibility of no observed resonance below the current lower bound on m* 
(~4 TeV for KK), one should be able to place a new lower bound on its mass as 
demonstrated in pp| . 

• Finally, the unique behavior of the invariant mass distribution of the KK model below 
the resonance may provide hints for the existence of such a resonance even if it is 
beyond the LHC reach. 

The preliminary analysis presented here is based on results from the combination of 
PythiaS and the new MoSES software. The specific Kaluza-Klein model is now imple- 
mented in Moses. As declared, this is the first phase towards integration inside PythiaS 
as an internal process. 

Other plans include interfacing this generator to the ATLAS detector simulation soft- 
ware for continuing the systematic analysis, generalize and re-order this framework as can 
be expected from the next releases and, enhance this framework with more BSM processes. 
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Appendices 



A Derivation of the KK tower for a 5(i real massless 
scalar field 

To demonstrate in practice the origin of the so called KK tower and the mass relation to 
the size of the parallel ED, i?, we can start by thinking of a toy model where we have a 
massless scalar field in flat 5(i space-time. The metric tensor is qab = (1, —1, —1, —1, il), 
assuming a space-like ED where we denote the ED coordinate by z with x being the regular 
four-coordinate pO|,|31|. 



First, we assume that the field $(x, z) satisfies the 5d Klein-Gordon (KG) equation 

$(x, z) = {d^d" - dl) $(x, z) = Q (37) 
A second assumption, followed from the form of the KG equation, is to separate the 

oo 

dependencies (of the x and z coordinates) and take $(x, 2) = ^ (f)n{x)ipn{z) . The next 

n=0 
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assumption is that ipniz) satisfies the equation 

dlM^) = (38) 
We also assume the orthogonahty of the different z dependent parts ze, 

j dztprn{z)tpn{z) = 5mn (39) 

Zl 

and finally, we write the bd action of the field with the boundary conditions (BC) taken to 
be [i}k{z)d-,i^n{z)]ll = 

Z2 

S = j d^x j dz^ (40) 



Zl 



By integrating over the extra coordinate z, we reduce the 5d action to an effective Ad one 

oo „ 

^ = J2 ^'^2 l-^^'^-^''^- - "''.'Pi] (41) 

n=0 

This action can be simply understood as the sum of distinct actions of 4d scalar fields 
(j)n{x) with different masses labeled by the integer index n. In fact, the various masses we 
observe in 4(i correspond to quantized values of the momentum along the extra coordinate 
z since 

= p2 = q^bPaP^ = Pf^P^ - pI 



(42) 



where 7714^^ is the "observed" mass in Ad. The quantized momentum in the ED is Pz = idz 
so by using Eq ^ we see that 

"mf] = rUn (43) 

where it is more appropriate to identify also the mass in Ad with the integer index n = 
0,1,2,.... The last term is often called the KK tower of excitations. This is a striking 
result since we obtained a mass property in 4d by starting from a fiat, 5d space-time with 
a massless field. Yet, this doesn't explain the origin of the mass relation to the size of the 
extra dimension R. To obtain this relation, we have to look at the specific ED topology 
and the BC. Looking on Eq |3^, we can expand ipn{z) in exponential wave-eigenfunctions 

^^{z) = C+e''""^ + C-e-'"^"^ (44) 
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where are constants. If we postulate that the ED is compact such that it is curled into a 
one- dimensional sphere of radius R then the specific translation invariance z ^ z + 2ttR 
of ipn{z) implies that the BC are periodical ze, 

ipinR) - ij{-7rR) = ^ ^2n^Rm„ _i = q ^ = - (45) 

R 

Therefore, Eq ^ translates to mf'^ = ^ where we see that the massless mode does exist as 
required, ie, = 0. However, in Eq ^ we saw that the mass we observe in 4d consists of 
the KK excited term ^ but also from another fixed term mo which is identified as the Ad 
state. We know that this term can be non-zero and that it is the Higgs mechanism which 
is responsible for creating the Ad mass mo. Thus, we obtained the full expression for the 
KK tower in this specific topology. 

<j = ^<+{iy (46) 

If we also define a parity operation on the interval z G [— vri?, ttR] then we obtain the 
z — >■ —z mapping. This implies that there are two special, fixed points z = 0, and z = ttR 
which are left invariant by this Z2 operation when combined with the periodicity property. 
The eigenfunctions that build ipn{z) must now respect the discrete Z2 parity symmetry 
so the solution for ipn{z) given in Eq.(|^) has to be modified to either an even solution 
AnCOS (^-z) or and odd one i?„sin (^-z). The massless mode can belong now only to the 
even solution. This geometry is called the /Z2 orbifold. 

The scalar field derivation is of course merely a single example of some illustrative 
methodological steps and assumptions concluding the origin of the KK tower and the mass 
relation to the ED radius R. 
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B The Moses software 

B.l Download, configure and build 

The source code can be viewed at [http: / /projects.hepforge.org/moses| , or downloaded from 



http: / / www.hepforge.org/ downloads/ moses| for the latest version. To prepare for installa- 



tion, make sure to have LHAPDF, ROOT, HepPDT and PythiaS installed (PythiaS should 
be configured with — enabled-shared). Make sure that the scripts root-config and Ihapdf- 
config are in your PATH (located in the bin directories of these packages). Copy the full 
paths to HepPDT and PythiaS base directories, these paths will need to be specified during 
configuration. To configure and build MoSES, 

wget http : / /www . hepf orge . org/ archive/moses/moses-<x> . <y> . <z> . tgz 
tar -zxvf moses-<x> . <y> . <z> . tgz 
cd moses-<x> . <y> . <z> 

source configure. sh — heppdtpath=/path/to/heppdt 

— pythia8path=/ path/to/ pythiaS 

make 

To clean the installation of compiled objects, run make clean. To complete purge of 
transient files, and to unset all the exported variables, run make cleanall. To unset all 
the exported variables, source configure. sh — unsetal. On the next logons, after the 
initial configuration, one has to run source ./config.sh to set the necessary variables 
for runtime and for recompilations. 



B.2 Basic usage 

All the executables go in the bin/execs directory and they can be ran directly from 
there. To see the list of executables, run Is -Irt bin/execs. To run a certain executable 
(executable_name), run $MOSESSYS/bin/executable_name. 

There is a special helper script for handling with running and re-making the vari- 
ous executables. To see all the helper script options, run ./helper — help. To see 
the list of executables, run ./helper — Is. To run some basic examples, one can exe- 
cute ./helper — run=usageExample_BuildProcess or — run=usageExample_LHAPDF or 
— run=usageExample_HepPDT. 
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All the output is essentially in the form of, either flat files (.dat etc.), or .root files 
(trees or canvases) or .eps files. All the executables output go by default to, data/, 
where this is determined automatically by the variable, $MOSESDATA/ = $MOSESSYS/data/ 
= /f ull/path/to/data/. This can be overridden by reexporting this variable, 

mkdir / your /path/to/the /new/datadir 

export MOSESDATA=/your/path/to/the/new/datadir 

Possibly, include the last line in the end of the conf ig. sh script which has to be executed 
once on every logon after the initial installation. To see the data in its default location, 
run Is -Irt data/ or in general. Is -Irt $MOSESDATA/, where the later should be used 
in the case of a user defined location. 



B.3 Structure 

In general, the basic structure is very simple so it can be very easy to either run the ex- 
amples and to allow further personal development on top of the existing structure. The 
common area (under the moses-<x> . <y> . <z>/examples) contains common code so it can 
be linked against any package (also under the examples directory). Each package (one 
directory) should correspond to a single model. In other words, a single hard process with 
the specific related tasks to be executed (except for the tasks that can be thought of as 
common to several packages). The user can copy code from these two areas for his own 
needs, add new libraries or add new packages just like it is done in the examples direc- 
tory. For a detailed picture of the MoSES software structure, see Fig The compilation 
flow is organized from the lowest level (core) to the highest level (models), so when stand- 
ing in one location which contains a Makefile and executing "make", everything that 
is at higher levels will also compile if changed. The moses-<x> . <y> . <z>/make .mk and 
moses-<x> . <y> . <z>/examples/common/coimnonmake .mk files contain global definitions at 
different levels to avoid rewriting these definition in each of the other Makefiles. Thus, 
either one or two of theses files have to be included in every Makefile depending on the 
context. Note that there can be many tasks placed under each model directory. That in- 
cludes oriented plugins and/or main programs like MC-programs (pythiaS-mains, pythiaS- 
plugins), theoretical-analysis programs, data-analysis programs etc. Therefore, there are 



some places where one can find specific higher levels than the ones sketched in Fig In 
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particular there are analysis directories that encapsulate several analysis tasks, some of 
which presented in this paper. 
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Figure 15: The MoSES software structure after installation. The basic compilation flow 
is from the lowest level Makefile (in moses-<x> . <y> . <z>j to the middle level Make- 
file (in moses-<x> . <y> . <z>/examples/commonj and to the highest level Makefile(s) (in 
moses-<x> . <y> . <z>/examples/model_i ). 
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